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ABSTRACT. Bovine coronavirus (BCoV) is a pathogen related to 
enteric and respiratory diseases in cattle worldwide. Enteric (BECoV) 
strains of BCoV are predominant in South America, and genetic 
investigations have been conducted to identify its relationship with 
isolates of respiratory origin (BRCoV). In this study, we used a BRCoV 
strain (BR-UEL11) derived from an outbreak of respiratory disease in 
feedlot cattle in southern Brazil, and compared the partial sequence of 
the polymorphic region of Spike (which was detected and sequenced 
by two distinct reverse transcription-polymerase chain reactions) with 
those of other BCoV strains. The phylogenetic relationship of BR- 
UEL11 with Brazilian BCoV, which is associated with calf diarrhea 
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and winter dysentery (enteric, BECoV; respiratory, BRCoV), and 
classical reference prototypes was analyzed. The analysis showed that 
the BRCoV strains from Brazil clustered with a clade that was distinct 
from most isolates associated with calf diarrhea (BECoV) and ancestral 
prototype strains such as Mebus, Nebraska, and LY VB. Furthermore, 
the BRCoV strains from Brazil clustered with a clade that contained 
recent strains associated with winter dysentery, showing 98-99% 
nucleotide identity with those strains. These results suggested that the 
Brazilian BCoV evolved from being solely enteric to a dual enteric and 
respiratory tropic virus. 


Key words: Bovine coronavirus; Bovine respiratory disease; 
Spike glycoprotein; Dual tropism; Phylogeny 


INTRODUCTION 


Coronaviruses cause respiratory and enteric disease in humans, birds, and other 
mammals (Siddell et al., 1983, Woo et al., 2009). The bovine coronavirus (BCoV) is a 
global pathogen that is recognized as the causative agent of enteric disease; however, the 
contribution of this pathogen towards the development of bovine respiratory disease (BRD) 
was underestimated compared to those of viruses that are primary agents of BRD, such as the 
bovine herpesvirus 1, bovine viral diarrhea virus, bovine parainfluenza-3 virus, and bovine 
respiratory syncytial virus (Fulton et al., 2011). However, recent reports have indicated that 
BCoV might also be a major pathogen of BRD (Lathrop et al., 2000; Boileau and Kapil, 2010; 
Hick et al., 2012). Moreover, isolates of BCoV associated with neonatal and adult diarrhea 
are referred to as BECoV (bovine enteric or enteropathogenic coronaviruses), whereas strains 
isolated from nasal secretions and lung of cattle affected by pneumonia have been classified as 
bovine respiratory coronaviruses, BRCoV (Lin et al., 2000; Boileau and Kapil, 2010). 

BCoV is a ubiquitous betacoronavirus within the subfamily Coronavirinae, family 
Coronaviridae (ICTV: http://www.ictvonline.org/virustaxonomy.asp). The virus is enveloped 
and pleomorphic, with a single-stranded positive-sense RNA genome. It has two major outer 
glycoproteins, namely, Spike (S), which is further divided into domains S1 and S2 (Cavanagh, 
1983), and hemagglutinin-esterase (HE) (Hasoksuz et al., 2002). S1 is the receptor-binding 
domain responsible for viral attachment and entry into the host cell, whereas S2 is related 
to membrane fusion. Therefore, S glycoproteins and HE are major neutralizing antigens of 
BCoV (Saif, 1993; Lin et al., 2000). Mutations in the S protein, mostly in the S1 domain, 
are associated with variations in the host range, tissue tropism, pathogenicity, and antigenic 
differences among strains (Saif, 2010). 

Worldwide, BCoV is associated with BRD in cattle of all ages, shipping fever, calf 
diarrhea, winter dysentery (Boileau and Kapil, 2010), and dysentery during the warmer 
seasons (Park et al., 2006; Decaro et al., 2008). Consequently, BCoV can simultaneously 
cause enteric and respiratory infections in the same animal (Park et al., 2007), although genetic 
markers for distinguishing between strains with both types of tropism were not identified 
during molecular analysis of the S glycoprotein (Kanno et al., 2007). Phylogenetic analyses 
based on the S-encoding gene sequence was performed in South Korea (Park et al., 2006), 
Japan (Kanno et al., 2007), and the USA (Fulton et al., 2013), which revealed that ancestral 
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reference strains of BCoV that are exclusively associated with enteric manifestations clustered 
within a distinct phylogenetic group compared to the recent isolates associated with enteric 
and respiratory diseases. 

Since the first description of a Brazilian BCoV associated with calf diarrhea, most 
Brazilian strains of enteric origin have now been described and characterized by molecular 
methods (Jerez et al., 2005; Takiuchi et al., 2008; Stipp et al., 2009; Ribeiro et al., 2016); in 
addition, assessment of major public databases did not reveal any description of BRCoV in 
other South American countries. This study investigates the recent molecular evolution of 
BCoV by addressing the phylogenetic relationship between the enteric and winter dysentery 
strains based on a polymorphic region of the Brazilian BRCoV. 


MATERIAL AND METHODS 


The BRCoV strain (BR-UEL11) was derived from nasopharyngeal swabs collected 
during an outbreak of BRD in feedlot cattle from the State of Parana, southern Brazil. Nucleic 
acid purification was performed as described (Boom et al., 1990), and the presence of BCoV 
was confirmed using a semi-nested polymerase chain reaction (PCR) targeting the N gene 
(Takiuchi et al., 2006). For genetic analysis, a reverse transcription (RT)-PCR assay using 
the SPK5 pair of primers (Takiuchi et al., 2008) was used to amplify the sequence encoding 
the nucleotides 1256 to 1904 based on the S gene of the Mebus strain (accession number - 
U00735.2), which includes the polymorphic region (nucleotides 1368 - 1776) (Rekik and 
Dea, 1994). All amplicons were visualized on an ethidium bromide-stained 2% agarose gel. 
The PCR products were purified using the GFXTM PCR DNA and gel band purification kit 
(GE Healthcare, Little Chalfont, UK) and sequenced (ABI3500 Genetic Analyzer, Applied 
Biosystems) using the forward and reverse primers with the BigDye™ Terminator v3.1 cycle 
sequencing kit (Applied Biosystems, Foster City, CA, USA). 

The phylogenetic analyses compared the BRCoV strain, BR-UEL11 (GenBank 
accession No. KP208779), with 67 Brazilian strains of BECoV and BRCoV, which were 
available in the GenBank. Furthermore, reference strains of respiratory origin, such as OK- 
0514-3 (AF058944) (Chouljenko et al., 1998), 232NS (DQ320763.1), the winter dysentery 
strain KW 10, and the classical calf diarrhea prototype enteric strains Mebus (U00735) (Stair et 
al., 1972), Nebraska (JQ741969), and the Brazilian LY VB strain (AY 606205) (Brandao et al., 
2006) were included in the study. To include all available data (from the Genbank) pertaining 
to the Brazilian BCoV, phylogenetic analysis was performed by evaluating the partial 
sequence of the polymorphic region (nucleotides 1388-1712). The phylogenetic relationship 
was constructed using the MEGA 7.0 software. The nucleotide sequences were aligned with 
Clustal W, and the phylogenetic analysis was performed using maximum likelihood’s Tamura 
3-parameter model and 1,000 bootstrap replications. 


RESULTS 


The phylogenetic analysis of nucleotide sequences from the polymorphic region 
revealed two groups that form the major clades | and 2 (Figure 1). Clade | consists of the 
BECoV strains from Brazil (USP- 2 to 14) derived from stool samples of calves affected with 
the winter dysentery WDBR-01 strain, and the ancestral calf diarrhea reference strains LY VB, 
Mebus, and Nebraska. A deletion of 18 nucleotides from positions 1577 to 1594 of the S 
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gene, which is exactly at the polymorphic region (Brandao et al., 2006), was present in all the 
Brazilian strains and the reference LY VB strain within clade 1. 
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Figure 1. Phylogenetic relationship, based on the partial sequence of the polymorphic region (nucleotides 1388- 
1712), of selected Brazilian strains of bovine coronavirus with prototype strains of BECoV and BRCoV. The 
nucleotide sequences were identified by GenBank accession number and strain name. BRCoV strains are denoted 
by a black diamond and BR-UEL11 with an arrow. The tree was generated using the maximum likelihood method 
with bootstrap values of 1000 replicates. 
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Most of the Brazilian BRCoV strains clustered in clade 2, including the respiratory 
strain BR-UELI1 identified during the outbreak of BRD (Figure 1), and other respiratory strains 
from Brazil. This clade also contained the respiratory prototypes OK-0514-3 and 232NS, as 
well as the winter dysentery KWD10 strain. Additionally, two subclades, 2A and 2B, were 
formed by viruses that originated from Brazil. Subclade 2B consisted of the respiratory BR- 
UELI1 strain, winter dysentery strains of enteric origin, and the BR-UEL 1, 2, and 3 strains 
obtained from diarrheic calves (Takiuchi et al., 2008). The subclade 2A consisted of several 
winter dysentery strains of enteric and respiratory origin and USP-1 from calf stool; however, 
no nucleotide deletion gaps, as observed in the USP strains of clade 1, were identified. 

Calculation of the genetic distance among strains revealed that the BR-UEL11 strain 
shared 99 to 98% nucleotide identity with similar strains from Brazil within clade 2, but only 
91 to 90% identity with strains from Brazil in clade 1, which harbored the 18 nucleotide gap 
in their sequence. Compared to the prototype strains, BR-UEL11 showed a nucleotide identity 
of 91% with LYVB, 96% with Mebus, 97% with Nebraska, 98% with KW10, and 99% with 
the OK-0514-3 and 232NS respiratory strains. 


DISCUSSION 


This study demonstrated the phylogenetic relationship among BCoV strains identified 
and analyzed in Brazil predominantly from 2003 to 2010. In addition, a sequence of the S1 
polymorphic region of a BRCoV strain identified in Brazil in 2012 was included in the study. 
Previous studies have suggested that genetic changes in the gene encoding the S protein of field 
strains may interfere with the humoral protection elicited by vaccination because vaccines are 
generally manufactured using the old prototype strains as antigens (Hussain et al., 1991; Lee 
et al., 2010; Kanno et al., 2013). Currently, commercial vaccines against respiratory pathogens 
of cattle containing BCoV-derived antigens are lacking in Brazil. Vaccines containing BCoV 
were developed to prevent calf diarrhea, being formulated only in association with other 
enteric pathogens. However, empirical analysis based on molecular diagnosis (data not shown) 
suggests that respiratory strains of BCoV may be significantly prevalent in Brazilian feedlots. 

Older hypotheses posit that BCoV strains have genetically evolved from the ancestral 
enteric viruses to distinct respiratory and enteric strains, which are distant from the old 
prototypes responsible for calf diarrhea and winter dysentery (Hasoksuz et al., 2002). However, 
based on the analysis of the S gene of South Korean isolates, where strains associated with 
dysentery clustered with a BRCoV isolate and were distant from the classic Mebus prototype, 
it was hypothesized that the same strain could cause both respiratory and enteric disease, 
and that BCoV had evolved to possess dual tropism (enteric and respiratory) (Park et al., 
2006). To investigate this possibility, an experiment was performed with the South Korean 
winter dysentery strain KWD3, which showed that calves developed diarrhea without clinical 
manifestation of respiratory disease, whereas histopathology demonstrated mild to severe 
interstitial pneumonia caused by the KWD3 strain (Park et al., 2007). 

Phylogenetic analysis of the Brazilian strains revealed that the respiratory BR- 
UEL11 and WDBR-R-(x) strains clustered within a distinct clade that was distant from the 
old enteric BCoV prototypes. Clade | contained the old prototypes such as LYVB, Mebus, 
and Nebraska, whereas clade 2 contained more recent strains associated with winter dysentery 
and respiratory tract disease (Figure 1). Similar results were obtained upon analyzing the S 
gene of winter dysentery strains in Japan (Kanno et al., 2007), and with isolates of enteric and 
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respiratory origin from the USA (Fulton et al., 2013). These findings support the theory that 
this phylogenetic distribution resulted from the evolution of BCoV strains to possess dual 
enteric and respiratory tropism. 

In a phylogenetic analysis with strains from the USA (Fulton et al., 2013), one clade 
that contained more recent respiratory and winter dysentery strains further diverged into 
subclades, which varied antigenically. These differences were assessed by virus neutralization, 
using convalescent cattle sera raised against all the American subclades that cross-react with 
each other, and the Nebraska prototype strain from a distinct clade (Fulton et al., 2013). 
Thus, the observed antigenic differences suggest that vaccines manufactured with calf 
diarrhea prototypes, such as the Mebus strain, may have reduced efficiency against certain 
strains circulating in the USA. The same phenomenon might occur with certain strains of 
other countries such as Brazil. However, challenge studies in validated animal models are 
required to ascertain the correlation between strain-specific antigenic distances and efficiency 
of vaccine protection in cattle infected with the current Brazilian field strains of BCoV (Fulton 
et al., 2015). 

In conclusion, we have demonstrated that the respiratory and enteric winter dysentery 
strains of the Brazilian bovine coronavirus included in this analysis are distant from the 
ancestral LY VB, Mebus, Nebraska, and old field strains that clustered with prototypes related 
to calf diarrhea. These findings suggest that the BCoV strains of Brazil probably diverged 
from enteric to dual tropism, and currently the same strain is associated with both enteric and 
respiratory diseases. In addition, since the Spike protein is the major neutralizing antigen of 
coronaviruses, the differences in the S gene sequence between the Brazilian strains and the old 
BCoV prototypes highlight the importance of assessing the efficacy of vaccines raised against 
the current field strains in Brazil. 
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